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The investigations of phase transitions on the hier- 
archical lattices have begun with the development of 
the Migdal-Kadanoff renormalization group methodic in 
which such lattices emerged as approximants of the ordi- 
nary crystal ones 3 -. Later the variety of hierarchical lat- 
tices with noninteger dimensionalities were introduced 4,5 
thus being the models of fractal structures. As fractals 
are ubiquitous in disordered media such as dilute crystals 
and porous materials^, the studies of phase transitions in 
them are of great interest. While the structures of hierar- 
chical lattices are not random they have broad distribu- 
tions of coordination numbers and characteristic lengths 
being in this respect the useful models of random sys- 
tems. Indeed the studies of the second-order transition 
in the Ising model on hierarchical lattices have shown 
the dependence of the critical indexes on the structural 
characteristics^ in close analogy with that found in ex- 
periments on disordered crystals under the variations of 
disorder strength 9 . 

The influence of structural inhomogeneity on the first- 
order transitions in hierarchical fractals also has qual- 
itative similarity to that found in the numerical stud- 
ies of such transitions in the Potts models in dilute 
crystala 10 i 1:L i 12 i 13 and porous media 14 i 15 . Thus the ef- 
fect of transformation of the first - ord er transitions into 
the second-order ones found in Refs. Il(ll lll'Jl.ll lllol 
exists also in the g-state Potts model (q — 4, 10) on 
two types of hierarchical lattices with fractal dimensions 
d > 2 16 . This unusual phenomenon is hard to explain 
in the framework of the standard phenomenology, de- 
scribing the smearing of first-order transitions in random 
media (diminishing or complete vanishing of the jumps of 
thermodynamic parameters) as a result of appearance of 
inhomogeneous two-phase state near transition^. Such 
physical picture can not explain the nature of the insta- 
bility which occurs in the inhomogeneous systems and 
results in the divergence of correlation length and criti- 
cal susceptibility 10 i 11 i 12 i 13 i 14 i 15 . 

Meanwhile the experimental studies of transitions in 
liquid crystalsiS^ and antiferromagnet MnO&L, con- 
fined into porous media, corroborate the possibility of 
the transition order transformation under the influence 



of inhomogeneities. The transformation of the struc- 
tural first-order transition Oh — > D^h into the sym- 
metry forbidden second-order transition takes place in 
the magnetite Fej,Oi under Zn doping 21 . Also the 
change of the transition order is observed in mixed crys- 
tals (KBr)i- x (KCN) x 22 ^ at the ferroelastic transition 
from the cubic phase to the orthorhombic one. Such 
transition in the ideal crystals is always of the first- 
order type^i, but in (KBr)i- x (KCN) x it becomes of 
the second-order type at x = 0.65, 0.7— and x = 0.732 3 . 
with elastic module C44 going to zero at the transition 
poinljS 3 .. 

Thus the experiments and numerical studies of realistic 
models show that inhomogineity not only can smooth the 
first-order jumps of thermodynamic parameters but can 
also induce their second-order singularities. The elucida- 
tion of the nature of this phenomenon and development 
of its theoretical description could be very important for 
the understanding the mechanisms of phase transitions in 
random media as well as for the numerous practical ap- 
plications of disordered materials. The task posed before 
theory is to determine the classes of first-order transitions 
and types of random media in which the second-order sin- 
gularities appear and those in which usual jump smearing 
take place owing to the existence of intermediate inho- 
mogeneous phas o 17 i 25 . Another task is to determine the 
critical indexes for inhomogeneity-induced second-order 
transitions and to study their dependence on the struc- 
tural characteristics of random media. 

In solving these problems essential help can be ob- 
tained from the studies of spin models on such simplified 
imitations of real inhomogeneous systems as hierarchical 
lattices, the thermodynamics of which allows exact ana- 
lytical description in some casesi&. Indeed the results of 
Ref. along with the investigations of Potts model on 
random scale-free graphs 2 ^* 2 ^ are the only analytical ev- 
idences of the possibility of transformation of first-order 
transitions into second-order ones in inhomogeneous sys- 
tems. 

Here we may note that the ch oice of Potts model 
for these studies in Refs. lilllTlTIT2ll3ll4ll5lltJ26l27l 

as the simplest model having first-order transitions on 
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FIG. 1: Graphical representation of the hierarchical lattice 
construction, m = 3 and n = 4. 



translationally-invariant lattices with d = 2, 3 is also stip- 
ulated by the existence of its numerous physical realiza- 
tions. Among them there are structural transitions in 
monolayers (q = 3,4), transitions in the cubic ferromag- 
nets in external field and in liquid mixtures (q = 3), see 
review Ref. and references therein. There are also a 
number of ferroelastic transitions, described by the Potts 
model, e.g. transition Oh — > D4J1 in spinels of NiCriO/± 
and Fe^Oi type and superconductors Nb^Sn and V3S1 
(q = 3) or charge ordering transitions in Yb^As^ and 
alloys of Mg 3 Cd type (q = 4), see Table IV.4 in Ref. H 
Here we consider transitions in q - state Potts model 
with arbitrary q on the two-parametric family of hierar- 
chical lattices with fractal dimensions d > 1 and average 
coordination number z between 2 and 4. Using analyt- 
ical approach differing from that of Ref. 

El 

it is possi- 
ble to show that second order transitions take place at 
all q, d and z, to obtain analytical expressions for criti- 
cal indexes, to study their dependence on lattice struc- 
tural parameters and to establish the scaling relations 
among them. The results obtained allow to suggest the 
structural criterion for the appearance of inhomogeneity- 
induced second-order transitions. Finally we discuss the 
application of present results to the description of criti- 
cal anomalies at phase transitions in dilute crystals and 
substances confined in porous media. 



I. GEOMETRICAL CHARACTERISTICS OF 
THE HIERARCHICAL LATTICES 

The procedure for construction of the hierarchical lat- 
tices we will consider here is depicted in FigQ] It consists 
in sequential substitutions of every bond by the n > 2 
chains having m > 2 bonds. 

At every step the number of bonds in the lattice be- 
comes B = mn times greater, so at fc-th step we have B k 
bonds in it. The number of sites Nk one can find from 
the recurrence relation 

N k = BNk-i -B-n + 2 

from which we have under the initial condition Nq = 2 



N, 



B-n 
B - 1 



-B h 



n-1 
B — 1 



1. 



(1) 



Introducing the average coordination number for infinite 
lattice, z = lim W— we get from Eq. 



z = 2 



B — 1 
B-n' 



(2) 



Thus z can vary from 2 (at m — > 00) up to 4 (at m = 
2,n — ► 00). At k-th step the largest distence between 

lattice cites is m k so Nk ~ (m fc ) d at k — > 00, where 

In 5 

d = 1 

mm 

is the fractal dimension of the lattice. Apparently, 1 < 
d < 00. 

Considering z and d as independent parameters from 
inequalities m > 2, n > 2 we get 



2 < z < 4(1 - 2~ d ) 



(3) 



Let us consider the distributions of cite coordination 
numbers in such lattices. The coordination numbers ac- 
quire the values Zk = 2n k and the number of sites having 
them in the lattice of the level I is 



s k = (B-n) B l -"- L ,0<k<l-l 

Also two basic sites have the coordination numbers zi — 

1 

n , so si = 2. One can easily verify that s& = Ni and 

fe=0 

/ 

^2 SkZk = 2Bi. Hence we have the following distribution 

fe=0 

function for coordination numbers in the infinite lattice 



W(z)= lim ^2^-5 (z-z k ) 

fc=0 



B - 1 
B 



2n k 



Thus W (z) has power law dependence with exponent 
greater than 1. Analogous distributions but with more 
dense sequence of Zk (z^ — k) are used for modeling the 
scale-free random graphs 26,27 .. Potts model with q > 1 on 
such graphs has second - order transition with mean-field 
singularities if exponent in W (z) is less than 3^&, which 
is related to the divergence of moment (z 2 ). In our case 
all moments {z r ) with r > d/ (d — 1) diverge, while at 
lower r 



(z r ) = 2' 



B-l 

B-n r 



Yet we show below that on the hierarchical lattices con- 
sidered the transition is always of the second-order type 
and the properties of W (z) influence only the values of 
the critical indexes. 
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II. RECURRENCE RELATIONS, 
THERMODYNAMIC POTENTIAL AND 
CORRELATIONS IN POTTS MODEL 

The partition function of the g-state Potts model on 
the hierarchical lattices considered can be obtain intro- 
ducing for each bond the factor 



a>i,l+i — 



a'21 + (g ~ 2) cj 
9-1 



, n n 



+ (q - 2) 



C/+1 



a 2l 




q - 




a 2l 




q 


- 1 



(9) 



According to Eq. (0J the initial conditions for these 
relations arc 



Z (a, a') = exp 



K5a,a' + 2 (<W + <W) 



(4) 



ai,o = e 



a 2 ,o = e K + g-2, 



6 = e 



/t/2 



co = e K - 1 



(10) 



Here a, a' are the Potts spins on the connected cites, 
K = J/T, h - external field. Summing the obtained 
expression over the values of Potts spins a = {1, 2, q} 
we get the full partition function. 

Introducing the partial partition functions, Zi (a, a'), 
for the level I lattice summed over all spins except cr, a' 
on the basic sites, one can get the recurrence relations 
for the m 1 ^? 



Solving Eqs. iJSJ one can find Potts partition func- 
tions with different boundary conditions on the basic 
lattice sites. Thus we get for free boundary conditions 
adding the missing fields h/2 to basic sites 



Z l+1 (cr, a') = 



Zi 



exp - (1 - n) (S a s + cv,i) 



(7,(7' 

e h a u + (?-!) (a 2 i + 2e' l / 2 6 / ) 



(11) 



( 5 ) 

Thus the next level partial partition function is ob- 
tained as m-th power of the matrix Z\ (cr, cr') with sub- 
sequent exponentiation of each its element to the n-th 
power. The exponential factor in Eq. (JSJ deletes ex- 
tra powers of exph(5 cr ^ + 8 a ^i) /2. Using Eqs. I@] [SJ 
one can find full partition function in the thermodynamic 
limit I — > 00. 

From Eqs. I@] EJ it follows that Zi (cr, a') can be rep- 
resented as 

Zi (cr, cr') = aiicv.icv,! + a 2 i (1 - S^i) (1 - cv,i) /(q - 1) 

+Q [S^ ~ 5as5a>S - (1 - <W) (1 - /(<? - 1)] (6) 

+h (8 a s + CV,1 + CVlC^i) . 

From Eq. EJit follows that the matrix Zi (cr, cr') has two 
nondegenerate eigenvalues, 



For the periodic boundary conditions identifying the ba- 
sic sites we have 



Z, 



(p) 



Zi ((7, cr) = an + a 2l + (q — 2) C\. 



(12) 



At last for the boundary conditions with fixed spins cr = 1 
on the basic sites 



.(i) _ 



= Zi (1, 1) = an 



(13) 



A±i = \ (ai,i + oa,i) ± Y j (ai,i - a 2 ,;) 2 + («-!) & f (<0 

and the eigenvalue c/ with the degeneracy q — 2. The 
transformation of the coefficients in the Eq. [S]when the 
matrix Z\ (cr, cr') is exponentiated to the m-th power can 
be represented using Eq. 0in the form 



In the absence of long-range order all these partition 
functions give in thermodynamic limit the same values 
of the thermodynamic potential per site. At h ^ their 
determination is rather hard task, yet the potential and 
its field derivatives can be more easily obtained at h = 
by the analytical methods. This is sufficient for determi- 
nation of the transition order and description of critical 
anomalies. 

In this paragraph we obtain the expression for the ther- 
modynamic potential in zero field. In this case we have 
from Eqs. JEIEJ 

an =bi + ci a% = an + (q-2) q, 

so in zero field there are only two independent coeffi- 
cients. Introducing 



+ = »iA'!\ + A m ; ) + |K/-«2/), 
4 = 5 (AT/ + A m ; ) - f (a u a 2l ) , 



K, = In 



an 
bi ' 



(8) 



b[ = cih, 4 = cf . 

= (A^-A^)/(A +i -A_0. 

Then the recurrence relations for coefficients correspond- 
ing to Eq. (JjJJ are 



we get the well-known relation^ 

e Kl+1 = f (e Kl ) 
q(x-ir 



(14) 



/(*) 



1 



(x - 1 + q)" 1 - (x - iy 
The second recurrence relation at h = is: 



-h{n-l)/2 (pl\n 



h+i = gi+ibf, 



= #1-1 



91 



1) J 



oKi/r. 



(15) 
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From Eq. Ill5[l and the definition of A; we get the ex- 
pression for the thermodynamic potential in zero field: 

_ oo 

F=-T lim Nr 1 InZP = --T V B~ k \ng k . (16) 



k=l 



To derive Eq. I|16[l we take into account that KiB~ l — » 
when Z — > oo. 

The relations in Eq. i|15|) have one stationary point 
K = K C , 



(17) 



corresponding to the transition point. At A > K c Ki — » 
oo, and at X < A c A; — > 0. When |A" - A" c | < A c , then 
for sufficiently small I Ki varies slowly, 

e K t -K c ~i + K i( K _ Kc ) i (i 8 ) 

« = /'(e*-) = 

B ~ n? tyt^ tv- < B ' (19) 

(e^ - 1) (e^ + g - 1) K ' 

until it becomes much larger or much less than A c . The 
condition for validity of Eq. (|18fl is 



const 



(20) 



Here the constant is determined from Eq. (|18|) by the 
condition K\ c ~ A c at A < A" c and by the condition 
A; c > A c at A > A c . At Z > / c it follows from Eq. fT5|) 

i-i c 

A < A c , A; « gn- 1 /'" 1 - 1 ) (A^nVC— (21) 

A > A c , exp A, sa m"^"" 1 ) (m""/^" 1 ) exp A ic ) 

Using Eqs. (|18l) - (|21|l one can show that near transi- 
tion the potential 116[) has singular part proportional to 
B~ lc ~ | A — A c | 2 ~" with the heat capacity index 



a = 2 - 



m.B 
InK 



(22) 



Thus at A > A c and | A — A c | ~ K c assuming in Eq. 
(EJ) 



9k 

we get 



g c ^{e K c-l) B /(e K °/"-l) n ,Kl c 
9oo =(q m -V n ) n ,l>l c 



2(5-1) 



-T [ flc + ( 5oo - ffc )] . 



At /i = it is also easy to find the correlation function 
for the basic spins a, a' 

G = (<5 CT ,l<5<r',l) - (S a l) (6a', i) ■ 



On the Z-th level lattice we have 



G, = 



Zi (1,1) 



g-1 e Ki -l 
g 2 e Kl + g - 1 



At A < A c , Z > Z c it follows from Eq. J2TJ 

G, « A ; /g 3 « g- 2 n-V(™-D (jf lo n 1 /( m - 1 Vg) m,_, ° 

~ exp (— m'/£) 

where the correlation length 

lnm 



£ ~ m<° ~ (A c - A)"" , 



In; 



(23) 



The correlation length index v obeys the scaling relation 

dv = 2 — a 

At A > A c from the second relation in Eq. (|21|l it follows 
that Gi goes to a constant at large I, 

Gi » (g - 1) /g 2 + A exp (-const ■ n l ~ la ) , 

but in this case the characteristic length of Gi varia- 
tions is also proportional to m lc , as n l ~ lc = (m l /m' c ) 
Thus on the lattices considered the Potts model has al- 
ways power low singularities of the heat capacity and the 
correlation length, with exponents obeying the usual scal- 
ing relation. In the next section we consider the critical 
behavior of the order parameter and the critical suscep- 
tibility. 



III. ORDER PARAMETER AND CRITICAL 
SUSCEPTIBILITY 



The spontaneous order parameter of the Potts model 



is 



lim 

; — >oo 



Ni-l 

QNf 1 £ (6 aiA )-l 

i=l 

g-i 



(24) 



Here (6 a ,i) means the average with boundary condition 
cr=l on the basic sites, that is 



Ni-l 



1 dZ, 



(i) 



Z W dh 
**l 



7(1) 



h=0 



Z, 



(1) 



(25) 



Just the use of such symmetry breaking boundary condi- 
tions allows to obtain a nonzero n at A > A c and h = 0. 
Indeed in zero field and for free and periodic boundary 
conditions (6 ai ,i) = 1/g due to the permutation symme- 
try of a values. This results, particularly, in the relations 

Z<« = N^P/q = (Nt - 1) Z\ p) /q 
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Inserting here Eqs (|llfl . (|12fl . we get 

A +i = {Ni - 1) A +; /g (26) 
A_, + (q - 2) Ci = (Ni -l)(q- 1) X-i/q 

Here A±z are the field derivatives of the eigenvalues in 
Eqs. at h = 

A±; = - (di; + d 2i )±^ (2 - (?) (di; - d 2 ;) + 2 (<? - 1) 6; 
From Eqs. JTSJ) (EH) - © we have 



/i = lim 



gJVfV+i+?-2 



i^oo 2 (g — 1) 
V+i = [ou - a 2 / - {q - 2) cj] /an. 



(27) 



Let us introduce one more linear combination of the 
derivatives linearly independent with (p + i and those in 
Eqs. HU), 

f-l = {an - a 2 i + qci) /an 

Then, differentiating Eq. © over h and using Eq. 1(26(1 
we get the reccurence relations for the vector ipi — 



where 



tpi = Titpi-i + u h 

i _l q - 2 
u+i = 1 - n H 7T— 



(28) 



n(JV,_i-l) 



(g - 2) (m - e;) - 2 (g - 1) (m - 1) e -^/™l(29) 



u-i = + (JVj - 1) 1 - -i (30) 

e; \ ei J 

~_n f ei [2 + (g - 2) mi?,] e; (g - 2) (1 - m^) \ . 
' q'\ 2e{(l-m*0 ej (g - 2 + 2mi9,) / dA; 

e ; = exp [Ki-i - Ki) , e[ = exp (ifj-i - K t ) , 

z? ; = [exp (Ki/n) - 1] / [exp (K^) - 1](32) 

The solution of Eq. is 

i-i 

<Pi = TiTj-i— fiy + y^Tjfi-i— ffc + iitfc +tt; (33) 
fc=i 

where (p = (1, 1). 

Let us consider the asymptotics of ipi at I — > oo near 
the transition. Here TJ can be approximated as 



lim TiA<L 



T„o = lim TiJ > L 



(34) 



Then 



fc=l 
i 



or in simpler form 



J- S^T, 



Here 



-B 



lim lim 



(35) 



lim it/_E> 

I— >oo 



<p c =(l-B 1 f c j U c ipoa 

From Eqs. JUJl-CSJl we get 

n ( e c [2 + (q - 2) toz? c ] e c (g - 2) (1 - mi? c ) 
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2(1- to?? c ) (g - 2 + 2mi? c ) 

e c = exp [if c (n — 1) /n] 
tf c = [exp (tf c /n) - 1] / [exp (K c ) - 1] 
2g-2 



36) 



z q 



It is important for the following that the eigenvalues 
of the matrix T c in Eq. 1)36(1 are real and are less than B 
all m > 2, n > 2, g > 0. Indeed, the largest eigenvalue 
of f r is 



2A max /n = to + £ - p + \J {m + e - pf - 4me (38) 

q-2 



p=-(l-e) 

q 



(to — 1) 



2 

£ = e c i? c < 1 



Using Eq. (|17|l for iT c one can show that < p < 
(to — 1) (1 — e) at all to > 2, n > 2, q > and this 
guarantees the reality of T c eigenvalues and the inequal- 
ity A 

max 

< B. 

The expressions for and Woo differ in ordered and 
disordered phases. Thus at K < K c it follows from Eqs. 
(ED, ED, (El, El 



n f 2 q-2 
°- q\2 q-2 

At K > K r 



q-2 m- 1 2 / 1 
g to z \ 1 











o 


J ' 


"oo = 





so the part in Eq. 1(35(1 containing Uqo (Voo) is zero and 
we have 

V+l = B l - l °e+f l c ° {<p - <fc) + B l <p c+ 
w iVi [consi (A max /B)' c - (g - 2) /g 

Thus we have for the order parameter in Eq. ((27|) 

M - (A max /S) /c ~ {K - K c f , 



k=l c + l 



p = MB/x max ) >Q 

In k 

(39) 
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Let us now consider the critical susceptibility at h = 



X = lim N- 1 



Z 



Z, 



(/)' 



z 



= lim (JVj-V+l - iV ; / g 2 ) + 2/g (40) 

t— *oo 

Here we used the relations in Eq. (|26fl and introduced 



a>u + (q - 1) (a 2 i + 26^ /<?A +; . 



Defining 

-0-; = 



(g - 1) I oii - 2bt) + a 2 i +q(q-2)ci / q\-i 



and differentiating Eq. twice over h we get the recur- 
rence relations for the vector tpi = (ip+i,ip^i) 



where 



Pi=B 



Xl (1-10/(9-1) 
(q-l)yi l-yi 

e Ki(n~l)/n _|_ _ ^ 



(41) 



(42) 



Xl = f e ^/ n + g -i 



yi 



1 



e A "< + q - 1 
e if, («-!)/« _ i 

e K > - 1 ' 



and at 2f c | ~ K c and large Z has the form 
B - 1 Nf ( 1 



B g 2 



cA 



2/ 



0(AV 



(43) 



The matrices has eigenvalues B and B (xi — yi) < B. 
The contribution to proportional to Nf is the right 
eigenvector for all Pi corresponding to the eigenvalue B. 
Using the last circumstance and the approximation for 
Pi, analogous to that in Eq. I|34|l . we get from Eqs. H4U|) 
- Eq. (|43[1 near transition 

X-(\LJB) 1 ^\K-K c \-\ 

7= (21nA max -lnB)/ln K . (44) 

Apparently the usual scaling relation is valid, 

a + 2/3 + 7 = 2. 

Thus the Potts model on the lattices considered always 
undergoes a second-order phase transition with power- 
law singularities of thermodynamic parameters. 

IV. CRITICAL INDEXES 

Here we consider the dependence of critical indexes on 
the lattice parameters and number of Potts states q. The 



inequalities k < B, A max < B (see Eqs. (ftTHl . an d 
scaling relations results in the following boundaries for 
the values of critical indexes in Eqs. it25|) . |J55J, 

u> 1/d, a < 1, /? > 0, 7 > 1. 

When g — > oo the indexes go to the limiting values in 
these inequalities. Indeed for q ^ 1 it follows from Eqs. 
03. <El © « (2/z)ln 9 , k « B, A rotKC » Be, so 



a 



1, f3^q- 2/in hiB 7»1, 



Simple expressions for the indexes can be obtained for 



n — > oo (d — > oo), when /iT c 

Amaa: ~ ^5 SO 



qn 



-l/(m-l) 



1, 



2-d, /3«1, 



d-2. 



Also at to — > co (d — > 1), when if c ss In to, k w n, 
p w gTO^r/6 , A max w B (1 - p/m), we get 



1 



d-l' 



d-2 

d~T' 



d 



d-l 



6m 2/(n-l) l n? 



Note that in these limiting cases only (3 depends on 
q and only when it is small. In general, considering the 
indexes as the functions of the observable fractal param- 
eters d, I, one can find that their values only slightly de- 
pend on the average coordination number in the interval 
2.5 < z < 4. Fig. J2J presents the indexes' dependencies 
on fractal dimension in the physically realizable interval 
1.5 < d < 3 at z = 2.5 and various q. At q = 1 the in- 
dexes describe the critical anomalies for the percolation 
of bonds randomly scattered over the hierarchical lattice 
with probability p — 1 — exp~ A , cf. Ref. 0HJ The equa- 
tion for the percolation threshold p c — 1— exp~ Kc follows 
from Eq. (fTTf) 



Pc = 1 - (1 



Let us note that all variety of the indexes of Potts 
model on fractals, illustrated in Fig. (J2J, has some defi- 
nite features. Thus a < for not very large q, the sus- 
ceptibility index is anomalously large (7 > 1.7), v shows 
the monotonous diminishing while (3 grows monotonously 
with a growth of the fractal dimension d. Note also a 
monotonous dependence of the indexes on q, the only ex- 
ception is j3. Probably, these properties of the indexes 
would be preserved in other spin models on the hierar- 
chical lattices at small enough d. 

We should also note that the indexes a and j3 at m = 2 
and some n and q were determined earlier in Ref. Il6i The 
expressions for a and /3 obtained here (Eqs. (|22|l and 
(13910 give the same values up to the calculation errors. 



7 




3,5 



3,0 



2,5 



2,0 



v 



A 



A 



~o-o-o-o-o 



_o_o-o- 



-A-A-A-A 



v-v-v-v-v-v-v 



V 



2,5 



1 ,5 



1,0 



0,5 



5 A 



-o-o-o— 



-A — A- 



■V—f, ~ 'A— A— A— A-^ 



3,0 



FIG. 2: The dependence of the critical indexes of the g-state Potts model on the hierarchical lattices with z — 2.5 on the 
fractal dimension d; □ — g = 1, Q~ Q = 2, A — g = 4, V — g = 8. 



ON A CRITERION OF TRANSITION 
ORDER TRANSFORMATION 



Considering hierarchical lattices as subclass of a wider 
set of inhomogeneous fractals which includes such real 
objects as porous materials and percolation clusters in di- 
lute crystals, one may attempt to draw some conclusions 
about their essential structural characteristics, which 
bring out the transformation of the first-order transitions 
into the second-order ones. Thus one can note that on 
the hierarchical fractals considered here this phenomenon 
does not depend on fractal dimension, at all 1 < d < oo 
the second-order transition appears. One can suppose 
that the feature that accounts for the transition order 
transformation in these lattices is their small average co- 
ordination number z < 4. Analogous situation appears 
in dilute 2d quadratic lattice, where dilution makes z < 4 
leaving the fractal dimension of the largest cluster d = 2 
up to the percolation threshold 6 . So, quite similarly, an 
arbitrary small concentration of vacancies transforms a 
first-order transition into a second-order one in this case, 



cf. Refs. At the same time, in the model of 

random media with z > 4 a first-order transition can be 
smeared but does not become a second-order one2^. 

It seems rather probable that the condition z < 4 could 
be a criterion of the transition order transformation in a 
rather large class of inhomogeneous systems with short- 
range interactions. Particularly, in dilute models on the 
simple cubic lattice (z — 6) with vacancies' concentration 
1 — x one may roughly estimate the average coordination 
number in the percolation cluster as z = 6x. Then from 
the condition z < 4 it follows that a second-order tran- 
sition would appear at x < 2/3. This estimate are close 
to the result of numerical studies, x < 0.7 ± 0.0 5 12 i 13 . 

Qualitatively the existence of the threshold coordi- 
nation number for the transition order transformation 
could be explained by the absence of the macroscopic 
regions with the largest coordination number in the sys- 
tems with low connectivity (low z). This makes ener- 
getically unfavorable the appearance of the macroscopic 
regions of the ordered phase near the transition points. 
In this case instead of creation of ordered regions only 
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spreading of the order parameter correlations is possible, 
which is characteristic for the second-order transitions 
and causes the singularities of the thermodynamic pa- 
rameters. Quite contrarily, in the structures with large z 
the phenomenological scenario of Ref. ^3 could be real- 
ized, which supposes the appearance and growth in vol- 
umes and numbers of macroscopic ordered regions near 
transition. Then both in the ordered regions and in the 
disordered ones the order parameter correlation length 
stays finite and the second-order singularities are absent. 
Yet in this case the first-order jumps can completely van- 
ish at large enough disorder—. 

VI. DISCUSSION AND CONCLUSIONS 

The results obtained here can be used to describe 
(more or less quantitatively) the critical anomalies near 
phase transitions described by the Potts model in the 
percolation clusters of dilute crystals or in substances 
confined in porous media, which have fractal structures 
similar to those of the hierarchical lattices considered. 
This possibility seems rather reasonable for q < 10 as 
then a < in conformity with the exact condition for this 
index in random systems^. The analytical expressions 
for critical indexes obtained here make their comparison 
with experimental data quite operational. 

We may note that the specific feature of the criti- 
cal indexes obtained here, which are determined mainly 
by the fractal dimension and only weakly depend on z 
at 2.5 < z < 4, allows to explain their slight varia- 
tions with the change of defect concentration in dilute 
crystals2ii2*ii*i2i±ii. Indeed, the fractal dimension of per- 
colation cluster, where transition takes place in such crys- 
tals, is almost constant up to the percolation threshold 6 , 
so slight changes of critical indexes result only from their 
weak dependence on the average coordination number. 
At the same time, one can expect more essential varia- 
tions of the critical indexes when the substance under- 
going transition is confined in porous matrices as their 
fractal dimension can vary in a wide interval—. 

Unfortunately now the comparison of the indexes ob- 



tained with their experimental values is impossible due 
to the absence of the detailed data on the critical anoma- 
lies near the inhomogeneity-induced second-order transi- 
tions. Mainly this is the consequence of the fact that only 
the recent theoretical work^iSiiiiiSiijiii 1 ^ have estab- 
lished the predetermination of this phenomenon. Mean- 
while, of old the numerous experiments are known, which 
find the second-order transitions in crystals which should 
have the first-order ones according to the Landau the- 
ory of p hase transitions^. In the light of the Refs. 
ll0yillll2llLM4lfl5l results one can consider such experi- 
mental data not as curiosity, but as a result of the pres- 
ence in crystals of a considerable amount of defects and 
impurities. The characteristic examples are Nb^Sn and 
ViSi crystals^ in which the distinction of the transitions 
from the second-order ones were established only on the 
samples of sufficiently high quality^. 

In conclusion we may note that, very probably, the 
inhomogeneity-induced second-order transitions are not 
the specific feature of the Potts model studied in 
Refs^WWW 6 and in the present work. It is 
naturally to suppose that this phenomenon can ex- 
ist at all transitions to a subgroup of high-symmetry 
phase, which in ideal crystals are described by the Lan- 
dau potential with a cubic invariant of order parameter 
components^. At the same time, there are neither ex- 
perimental nor theoretical evidences of the appearance of 
this phenomenon when the relation group-subgroup be- 
tween the symmetries of disordered and ordered phases 
is absent, that is for so-called reconstructive first-order 
transitions^. Probably, in this case the appearance of 
the inhomogeneity-induced second-order singularities is 
impossible. 
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